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Abstract 

Here is presented the calculation of the dynamic electrical conductivity of fully ionized, strongly 
coupled plasmas as a function of the external electric field frequency a>. The calculations are 
based on the the formula for the energy-dependent collision frequency which is determined by 
means of the Green function theory methods, as a sum over the Matsubara frequencies. The do- 
main of extremely high electron density: 10 21 <n e < 10 24 cm 3 , and for the temperature varying 
from 10 kK to 1.000 kK was examined. The real and imaginary parts of the conductivity for 
every electron density are presented in the generalized Drude-like form as a two-parameter func- 
tion of the frequency u in the region < to < Q.5u p , where to p is the plasma frequency. A good 
agreement between the obtained results and the existing theoretical and computing simulation 
data is shown. 
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1. Introduction 

The defining of the dynamic high frequency (HF) electrical conductivity of the strongly cou- 
pled plasma is an actual and very important problem because such characteristics of plasma as 
the HF dielectric function, the coefficients of refractivity and reflectivity, and skin depth are 
expressed just in terms of the dynamic electrical conductivity. However, until recently the deter- 
mination of the dynamic plasma conductivity in the wide ranges of the electron densities (n e ) and 
the temperatures {T) was not taken as the priority task. Namely, in the most cases the mentioned 
problem was not investigated as the main one oJ-Ql], and only few authors presented the results 
of the direct definition of HF electrical conductivity related to some selected values of n e 

and T. 

Because of that, a few years ago a certain research of the HF conductivity of the plasmas 
was started. It was thought that the previous researches of the transport properties of the dense 
strongly coupled plasmas free of any external fields and the plasmas in the presence of 

the constant magnetic field II 1 311 . gave some results which were useful in the case of the plasma in 
the external time depending fields. From this reason the research of the HF plasma conductivity 
was also devoted to the dense strongly coupled plasmas. The main aim was to examine the HF 
conductivity in the whole area of such plasmas, including here the region of the extremely large 
degrees of its non-ideality. This was especially important, since the existing literature data about 
the HF plasma conductivity lIB- Holl are just related to the extremely strongly coupled plasmas 
(10 21 cnr 3 < n e < 10 24 cnr 3 ). 
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The first results of that research were presented in ffl . where the new method of the cal- 
culations of the strongly coupled plasma HF conductivity is described. This method was tested 
on the plasmas with 10 17 cirT 3 < n e < 10 19 cirT 3 and 10 4 K < T < 5 ■ 10 4 K, in the microwave 
and far-infrared regions of a>, were <±> is the frequency of the external electric field. The devel- 
oped method uses the energy-dependent electron collision frequency which is determined as a 
sum over the Matsubara frequencies, using the Green function theory methods and the corre- 
sponding self-consistent calculations procedure, which is described in details in ifTll Il3ll . The 
final results, i.e. the HF conductivity <j{oS) = cr(u>; n e , T), as well as its real and imaginary parts 
o"«e,/m(w) = o-R e j m (aj; n e , T), were presented in the adequate, three-parameter Drude-like form. 
As the second step, by 01511 the same method with its improved numerical procedure it was ap- 
plied to the denser plasmas, with 10 19 cm 3 < n e < 10 2I cm 3 and 2 • 10 4 K < T < 10 6 K, in 
infrared, visible and near UV regions of co. 

Since, in lfl4l[l5ll the domain of low and moderately non-ideal plasmas was covered, the first 
aim of this Letter was to apply the developed method by the adequate numerical procedure, in 
order to determine of the HF conductivities of the strongly and extremely coupled plasmas, i.e. 
for 10 21 cm 3 < n e < 10 24 cm . This way would not only laboratory strongly coupled plasmas 
but also some astrophysical plasmas (for examples, the photospheres of some white dwarfs lHH]) 
be considered. 

Another aim is to compare the obtained results with the existing theoretical results and 
computer simulation data Bq- llOll . which refer to n e = 3.8 • 10 21 cirT 3 , 2.52 ■ 10 22 crrT 3 and 
1.61 • 10 24 cirT 3 . Such a comparison is very important because there is no data concerning the 
direct experimental measurements of HF conductivity of the dense strongly coupled plasmas [4], 
& 

Similarly to our previous papers jl4 Il5ll . the HF conductivity cr(w), as well as o~R e j m (u>), 
are presented here in the three-parameters Drude-like form, where one of these parameters is 
the static conductivity <x of the considered plasma, which is suitable for further consideration. 



However, while in 1114 11511 the values of cr could be taken from [13], the corresponding values 



for the region 10 21 cirT 3 < n e < 10 24 cirT 3 were especially calculated here. 

The calculations of o- Re (u>) and o- Im (a>) in all considered cases, for all taken values of n e and 
T, are performed here in the frequency range < a> < 0.5 -cOn, where a> p is the plasma electron 
frequency. The results are compared with the data from Il6l 4l0ll . 



2. Model 



2.1. Basic equations 

In this Letter it is considered a completely ionized plasma in a homogenous and monochro- 
matic external electric field ~t = Eo ■ exp{—iu>t}. According to lll4ll . the dynamic electric conduc- 
tivity of a strongly coupled plasma cr(w) = <j Re {(jj) + io- Im (aj) is presented by the expressions 
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where e, m and E are the charge, mass and energy of the free electron, t(E) - the relaxation 
time presented later in the text, p(E) - one-electron states density in the energy space, w(E) = 
[exp(J3E - flfi) + l] -1 - the Fermi-Dirac distribution function, p. - the chemical potential of the 
ideal gas of the free electrons with the density n e and the temperature T, and [5 - (ksT)' 1 . 

Comparing Eq. (Q]) by 10 — with the exp ression for the static conductivity <x of the fully 
ionized strongly coupled plasma from 111 lL 1 1 311 - one can see that 



cr(w = 0) = ctq 
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Consequently, from Eqs.(f2| and (O it follows the validity of the relations cr Re (u) = 0) = cr , 
(Ti m ((x) = 0) = 0, which have the sense of the conditions in the point o) = during the numerical 
calculations of cr Se (o») and cr Im {u)). The right side of Eq.© has the Lorentz-like form. Because 



of that the quantity t(E) was interpreted in 111 lL 1 1 311 as the corresponding electron relaxation time 



(total electron-ion plus electron-electron), characterizing the considered plasma in the absence of 
the external fields. Since the same quantity t(E) is used in the expressions ([T}[3]) for cr(a>) just the 
good understanding of the static conductivity <x as a function of n e and T has a very important 
role in defining the strongly coupled plasma dynamic conductivity. Because of this fact our 



research was started in (see 114IH in the region of the plasma electron density n e < 10 19 cm 3 , 



where the basic method of the determination of the <x from 
III- 



1 1 L 1 1 3tl was experimentally verified 



Accordingly to 11 111 1311 . the relaxation time t(E) is defined in the self-consistent approxima- 



tion by the following relations: 



T = 



\nmn e e' 



\lHmE/h 



a,v 



(q) ' 



(5) 



where a is the label for the plasma species, Z fl - the species charge number, Tl av (q) - correspond- 
ing polarization operator in the random-phase approximation (a simple loop), s v - the dielectric 
function expressed in terms of these quantities, and the v summation is extended over the Mat- 
subara even frequencies . Here, the polarization operator II,Q-) y (^) of non-degenerate j'-th ion 
component is equal to pn^dyQ, and the RPA polarization operator Tl ev (q) of partially degenerate 
free electron component, namely 



m ( kf\ 
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where k F = (3w 2 n f )^ 3 , Q v = 4nvm/h 2 f3q, v = 0, ±1, and 6 = yzm/p> ■ \hkf) 

The main numerical problems, which appear during the calculation of the static conductivity 
ctq, are connected to the defining of the relaxation time r. Namely, the number of the members in 
the y sum in Eq. (0), which has to be calculated individually, rapidly increases with the increase 
of the electron density n e . For example, at n e ~ 10 24 cirT 3 this number is order of magnitude 
~ 10 5 , while each member itself has to be calculated by a special procedure. It is clear that the 
determination of the HF conductivity o-(aj), requires the volume of the calculations several times 
larger than in the case of the determination of <t . 
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2.2. Parameterized form 

According to [141 and lfl5ll . the quantities o-R e and cr /m can be presented in the parameterized 
form, suitable for the comparison with the corresponding Drude-like formulas, namely 

1 (WTy)&2 

ctrJcj) = <t — — , <T /m (w) = (To — — , (7) 

\+(cjTlfk\ l+(cjT* Q ) 2 k 2 

where <x is given by Eq. ©, t* q is the effective relaxation time defined by <x = ^-t! and 
£1 = k\(to) and k 2 = k 2 (a>) are parameters, which describe the influence of the deviation of the 
dw/dE from the 5-function [13]. When the electron component of the considered plasma can be 
treated as a highly degenerate gas of the free electrons, both coefficients k\ and k 2 tend to 1, and 



equation O reduce to the Drude-Lorentz model 01911 . Using the plasma frequency cj p and the 
dimensionless parameter fo p , defined by relations a> p = {Ann e e 2 /m) 1 ^ 2 , fo p = (o p t1 = AitaolOp, 
the parameterized expressions for cr Re (w) and 0/ m (a>) convert into 

, , 1 , , (a>/oj p )fo p k 2 

o- Re {u) = cr ■ - — — — , o- Im (co) = 0-0 ■ - — — — , (8) 

1 + (oj/aj p ) 2 f 2 p k 2 1 + (aj/to p ) 2 f 2 p k 2 

Our previous numerical analysis have shown that k\ and k 2 can be cast in the following approxi- 
mate fitting form 

ajbj(oj/o) p ) 

kj = kjo -a r — — — — -, j = 1,2, (9) 

1 + ajbj(a)/a> p ) 

where the adjustment parameters k\o, a\,b\ and £20, a 2, ^2 are determined numerically from the 
condition of the best relation between the 'exact' and fitted values of the factors k\ and k 2 . Let 
us notify that by just using the expressions (IT)-© f° r o~R e {tS) and cr /m (<y) we get the possibility 
to confirm whether Drude-like model (when A4 2 ~ 1) has the physical sense for the considered 
n e , T and a>. 



3. Results and discussion 

In this paper it was defined the dynamic conductivity <t(u>) of completely ionized hydrogen- 
like plasmas, when in expression (0 only two species (electron and ion with Z = 1) are taken into 
account. Here we had on mind that the using of the expressions (Q) and ([8]) for the determination 
of o-R e (ai) and cr /m («) require the knowledge of the static conductivity <x . From this reason the 
values of <x are especially determined here in the considered domain 10 21 <n e < 10 24 cirT 3 and 
lOkK < T < lOOOkK. These values are presented in Tables Q] 

Then, according to [14] and [15], the coefficients k\ and k 2 in Eqs. (0 and ([8]) are taken in 



the form (O, where they are expressed in terms of the parameters kjo, aj and bj, j — 1,2. The 
values of these parameters were also computed in the same domain (10 21 < n e < 10 24 crrT 3 and 
10000K < T < 1000000K) and presented in Tableland Table [3] 

The conduct of o- Re (cj) and o"/ m (w), calculated directly by Eqs. (f2]l and (O for 10 21 < n e < 
10 24 cm and lOkK < T < lOOOkK., is displayed in tables and figures in the ranges < to < 0.5 • 
u> p . The figures[T][2]0and[6] demonstrate the regular behavior of o-R e (a>), i.e. the convergence to 
the corresponding values of (To(n e , T) when 10 — > 0, and the existence of the interval of variation 
of 10 where o-R e (u>) is practically constant. We observe the tendency of this interval to decrease 
when temperature T increases. Similarly, the figures[3][5]and[6l demonstrate a regular behavior of 
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cr /m (w), i.e. the convergence to zero when to — > 0, and the presence of a maximum in the interval 
< oj < Q.5oj p . The obtained results, together with the data from lSl, |@], J3l and ifioll . show that 
the position of this maximum comes to 0.5^, when n e grows in the region n e < 5 • 10 22 crrT 3 , 
and decrease from Q.5cj p when n e increases in the region n e > 5 ■ 10 22 cirT 3 . 

We compare the behavior of cr Re (a>) and cr/,„(w), determined in this Letter, with the behavior 
of the corresponding quantities determined in: |Q] and |01, for F = 0.5 and r s = 1 ; J3l and JUl, 
for T = 0.5 and r s = 4; H for T = 10 and r s = 1 ; for n e = 3.8 • 10 21 cm 3 andT = 3.3- 10 4 K. 
Here T = e 2 /(akT) and r s = a/ao are the well-known non-ideality and Brueckner parameters, 
where ao and a = [3/(4jTn e )] 1 ^ are the Bohr and the Wigner-Seitz radii. The results of 01, 
[8] and [ 10] were obtained by using the molecular dynamics (MD) simulation method, while the 
results of J6[] and |9[] are theoretical. The data from [6] are obtained by two analytical expressions 
for the collision frequency (the main one, obtained within the memory function formalism, and 
the approximate one found by the solution of the high-frequency hydrodynamic equation derived 
from the macroscopic equations of electron motion). The corresponding curves in the example 
T = 0.5 and r s = 1 are shown in Figs.Q]and[3j and in the case T = 0.5 and r s = 4 - in Figs.|4]and 

H 

The figures Q] and [3] show that in the example F = 0.5 and r s — 1 our values of <x fi( ,(w) and 
cr /m (w) are mostly in good relation with the results obtained in |0] by the MD simulation method, 
except the region a> < 0.1 ■ a> p . The same conclusion is valid for the theoretical results obtained 
in ((j] for r = 0.5 and r s — 1 by means of the mentioned approximate expression. Connected to 
the theoretical results from (fj and obtained by means of the main analytical expression in 
the cases F = 0.5 and F = 10, for the same r s = 1, a good agreement with our results exists in 
the whole region of a). Therefore it cold be expected since that main analytical expression has 
significant similarity with our expressions. 

From figure[3] concerning cr Im (oj), one can see that not only our results, but both results of |0] 
significantly disagree with the results of the MD simulation from [7] in the region to < 0.1 ■ a> p . 
We believe that final conclusions, respecting this fact, can be reached only after some additional 
MD simulations, since in this region the behavior of of the static conductivity ctq as a function of 
n e and T has the dominant role, and till now we have been considering our values of cr as fairly 
reliable. 

In the example when F = 0.5 and r s — 4 the figure [4] shows a good relation of our values of 
o-R e (aS) with the values obtained in J(J and Jst], since only in the region a> < 0.1 ■ a> p there exists 
some appreciable difference (the largest differences are less that 10 percents) between our results 
and those of [6]. The figure[5] concerning cr /m (a/), shows a qualitative agreement (the differences 
are less or close to 20 percents) of our results with the results of J6|] and [8J]. 



4. Conclusions 



The results obtained in this Letter, together with the ones presented in 1 14] and lfl5ll . make 



possible easy determination of cr Re (to) and cr Im (aj) in very wide ranges of the electron density 
and temperature, namely 10 17 < n e < 10 24 cirT 3 and 10 3 K < T < 10 6 K. Consequently, they 
can be used for the theoretical research, as well as for the interpretation of experimental data in 
the cases of high pressure discharge, capillary discharge, shock waves, etc., where the strongly 
coupled plasmas, including extremely dense plasmas, are created. Apart of that, cr Sc (w) and 
cr /m (w) can be useful in some cases of stellar plasmas (e.g., white dwarfs). It is important that 
the presented results allow the determination of other optical characteristics of strongly coupled 
plasmas: the permeability, the refractivity, etc. The advantage of the method, which was used in 
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our research until now, is that it does not contain any empirical parameters, and have got internal 
possibilities for its further improvement. 

Let us emphasize that the used procedure is similar to the procedure for the determination of 
the conductivity tensor in the case of a strongly coupled plasma in an external static magnetic 
field [13]. This fact opens the possibility for the development of the method of determination of 
the conductivity of strongly coupled plasmas in an external high frequency electromagnetic field 
in the presence of the static (or quasi-static) external magnetic field. 

Finally, the method used in this Letter does not assume in advance the quasi-Drude shape 
for cr(a>), but the way of presentation of the final expressions for cr Re and cr lm automatically 
demonstrates the character of the deviations from Drude's expression, at least in the case we 
consider. 

Certainly, since the numerical procedure for the determination of cr Re {<x>) and C7 m (a>) which is 
used here is quite complicated, it is desirable to develop some other methods which would allow 
an easier calculation of these quantities. Apart of that, in the lower part of the considered domain 
of the electron density (n e < 10 20 cm -3 ) it is necessary to develop the method of determination 
of the dynamic conductivity which would take into account the presence of the plasma neutral 
component. 
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Table 1: The computed static plasma conductivity <tq as a function of the electron density n e and temperature T, 

[lo'cnm) -1 ]. 



T[K] 








n e [cm J ] 












5 ■ 10 21 


10 22 


5 ■ 1(F 


10 iJ 


5 ■ 10 2 j 


10 24 


10000 


28.56 


78.23 


122.1 


388.0 


677.7 


2579 


4016 


15000 


33.90 


86.68 


135.7 


410.4 


701.0 


2698 


5175 


20000 


40.17 


93.34 


145.2 


430.8 


723.2 


2721 


4886 


30000 


54.19 


108.2 


160.8 


464.6 


762.6 


2767 


5037 


50000 


85.37 


145.7 


198.4 


511.8 


824.2 


2816 


5070 


100000 


174.2 


257.6 


320.1 


636.6 


948.7 


2920 


5070 


200000 


382.8 


516.8 


606.6 


985.9 


1305 


3200 


5228 


500000 


1165 


1461 


1639 


2278 


2726 


4779 


6669 


1000000 


2806 


3393 


3727 


4829 


5531 


8282 


10430 




Figure 1 : Plasma dynamic conductivity in [CGS] units, as a function of (oj/tOp): for T = 0.5 and r s = 1 , with the results 
from 01 and Ql, real part. 
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Table 2: Parameters kio (first row for each T), a\ and b\ (second and third rows, respectively) as functions of n e and T. 



T[K] 








n e [cm 3 ] 










10 21 


5 • 10 21 


id 12 


5 • 10 i2 




5 ■ m li 


10 24 



10000 



15000 



20000 



30000 



50000 



100000 



200000 



500000 



1000000 



1.36558 
0.91630 
0.80560 
1.56925 
1.12789 
1.06607 
1.70187 
1.22186 
1.40904 
1.85027 
1.30102 
2.22529 
1.98074 
1.36853 
4.09404 
2.09206 
1.45707 
9.00543 
2.16088 
1.53680 
19.0300 
2.21152 
1.59191 
57.9900 
2.26912 
1.65679 
138.900 



1.14298 
0.42670 
1.46304 
1.21020 
0.47940 
2.22226 
1.31099 
0.63100 
2.20497 
1.50720 
0.90070 
2.19804 
1.74157 
1.13885 
3.00424 
1.95358 
1.31413 
5.95428 
2.07644 
1.43764 
12.0400 
2.16940 
1.54523 
32.3600 
2.20725 
1.58655 
76.0200 



1.11647 
0.36160 
1.71715 
1.14482 
0.35200 
2.84164 
1.19312 
0.40760 
3.36852 
1.33432 
0.61770 
3.11968 
1.58508 
0.94640 
3.14523 
1.85948 
1.21608 
5.25356 
2.02196 
1.37555 
10.2000 
2.13967 
1.51208 
25.8200 
2.19266 
1.56883 
59.2600 



1.04503 
0.11570 
9.66853 
1.07335 
0.17550 
7.27586 
1.09268 
0.20970 
6.89132 
1.12245 
0.24750 
7.66151 
1.21653 
0.37520 
7.67889 
1.51913 
0.79270 
5.95785 
1.81532 
1.13930 
7.99801 
2.03771 
1.39149 
16.9600 
2.12778 
1.49589 
34.3800 



1.02252 
0.05374 
29.2500 
1.04339 
0.09703 
17.5200 
1.06185 
0.13380 
13.5600 

I. 08882 
0.18120 

II. 4900 

I. 13706 
0.24640 

II. 7600 
1.34490 
0.53610 
8.79760 
1.67551 
0.96530 
8.30448 
1.96910 
1.30924 
14.9600 
2.08724 
1.44904 
28.2200 



1.03775 
0.00832 
447.500 
1.00815 
0.01560 
244.600 
1.01567 
0.02620 
147.200 
1.02599 
0.04943 
80.3700 
1.05221 
0.09391 
44.5000 
1.10857 
0.17830 
29.6900 
1.28276 
0.41640 
19.0000 
1.70479 
0.97700 
14.5800 
1.93392 
1.26324 
20.7800 



1.23393 
0.00572 
975.700 
0.96590 
0.00541 
1032.00 
1.02906 
0.01223 
459.700 
1.01296 
0.02285 
252.400 
1.02846 
0.05013 
117.400 
1.06780 
0.11230 
56.2000 
1.16478 
0.24720 
34.4800 
1.53717 
0.75090 
18.2200 
1.82624 
1.12668 
20.2900 
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Table 3: Parameters &20 (first row for each T), «2 and bi (second and third rows, respectively) as functions of n e and T. 



T IK] 








n e [cm a ] 








10 21 


5 • lO^ 1 




5 • W l 


10 iJ 


5 • 10 2J 


10 24 


10000 


1.22578 


1.09158 


1.07647 


1.02936 


1.01445 


1.03622 


1.23283 




0.85360 


0.43490 


0.38030 


0.11700 


0.05381 


0.00821 


0.005572 




0.73820 


1.23034 


1.42388 


9.08148 


28.4700 


451.400 


1033 


15000 


1.36175 


1.12968 


1.09216 


1.04821 


1.02833 


1 .00524 


0.9649 




1.05294 


0.44090 


0.34460 


0.17850 


0.09769 


0.01544 


0.005332 




96230 


2 16439 


2 57804 


6.67217 


1 6 6700 

J. VJ ,\J / \J\J 


246 000 


1 066 


20000 


1.45309 


1.19154 


1.11970 


1.06089 


1.04061 


1.01077 


1.02666 




1.13563 


0.56690 


0.37520 


0.21280 


0.13520 


0.02613 


0.01221 




1.26874 


2.20112 


3.31066 


6.26462 


12.6900 


145.600 


457.6 


30000 


1.55543 


1.32016 


1.20668 


1.07872 


1.05821 


1.01668 


1.0084 




1.19391 


0.81940 


0.55200 


0.24180 


0.18220 


0.04959 


0.02282 




2.01639 


2.10622 


3.14718 


7.16940 


10.6200 


77.8800 


249.7 


50000 


1.64423 


1.48078 


1.37339 


1.13410 


1.08691 


1.03422 


1.01839 




1.23540 


1.04367 


0.86190 


0.33860 


0.23370 


0.09454 


0.05031 




3.76200 


2.76140 


2.98249 


7.80417 


11.4000 


42.0900 


113.6 


100000 


1.71861 


1.62588 


1.56188 


1.32945 


1.21494 


1.06960 


1.04423 




1.29643 


1.18484 


1.10384 


0.71430 


0.47780 


0.17160 


0.1116 




8.46413 


5.51045 


4.84527 


5.80599 


8.93956 


28.6500 


53.33 


200000 


1.76372 


1.70824 


1.67192 


1.53184 


1.43625 


1.17662 


1.10373 




1.36115 


1.28090 


1.23103 


1.03249 


0.87480 


0.37330 


0.2281 




17.9300 


11.3300 


9.55919 


7.45556 


7.86760 


19.4000 


34.65 


500000 


1.79739 


1.76921 


1.74994 


1.68242 


1.63638 


1.45691 


1.34383 




1.40905 


1.37086 


1.34419 


1.24460 


1.17609 


0.88420 


0.677 




53.6500 


30.0600 


24.1100 


15.9400 


14.0400 


13.9000 


17.87 


1000000 


1.82998 


1.79479 


1.78467 


1.74212 


1.71534 


1.61265 


1.53976 




1.45044 


1.40480 


1.39017 


1.33254 


1.29409 


1.13771 


1.01851 




130.200 


70.5400 


54.8300 


31.8700 


26.2700 


19.5000 


19.17 
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Figure 2: Plasma dynamic conductivity in [CGS] units, as a function of (ai/u> p ): for T = 10 and r s = 1 with the results 
from 1 9], real part. 
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The present work, r=0.5 




Figure 3: Same as in Fig. [TJbut for imaginary part. 
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Figure 6: Real and imaginary parts of the plasma dynamic conductivity for n e = 3.8 ■ 10 21 cirT 3 , T = 3.3 ■ 10 4 K, as a 
function of (w/u) p ) frequency ratio, for < oj/oj p < 1, together with result from 1 10]. The value to/a/ p = 0.5 especially 
mark to remind that within our research the region < wfat p < 0.5 is examined. 
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